Odd-frequency Superconductivity on a Quasi-one-dimensional Triangular Lattice in 

the Hubbard Model 
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In order to clarify whether the odd-frequency superconductivity can be realized or not, we study 
a quasi-one-dimensional triangular lattice in the Hubbard model using the random phase approx- 
imation (RPA) and the fluctuation exchange (FLEX) approximation. We find that odd-frequency 
spin-singlet p-wave pairing can be enhanced on a quasi-one-dimensional isosceles triangular lattice. 



I. INTRODUCTION 



It is known that two electrons form Cooper pairs with 
gap function Acr^a-2{'>'^n,k) in superconductors. In gen- 
eral, the gap function depends on the Matsubara fre- 
quency e„, the combination of spins CTi and CT2, and 



the momentum k. 
Dirac statistics, the 



In accordance with the Fermi- 
sign of the gap function is re- 



versed by the exchange of two electrons: Ao-jo-a ^s) = 
— AcTaCTi {—iSn, —k). Based on the symmetrical properties 
combinations of three dependences, symmetries of the 
gap functions can be classified into four groups; (i) even 
parity in frequency space (even- frequency) , spin-singlet, 
and even parity in momentum space (even-parity) labeled 
as ESE pairing state, (ii) even-frequency spin-triplet odd- 
parity (ETO) pairing state, (iii) odd-frequency spin- 
singlet odd-parity (OSO) pairing state, and (iv) odd- 
frequency spin-triplet even-parity (OTE) pairing state. 
Almost all of the superconductors including high-Tc 
cuprates belong to ESE pairing state. ETO pairing state 
is realized in some special superconductors like Sr2Ru04 
or UPta. In contrast to these even- frequency pairings, 
odd- frequency pairings are not familiar. 

In 1974, the possibility of realizing the odd-frequency 
pairing state was first proposed by Berezinskii in the 
context of ^He, where the odd-frequency spin-triplet hy- 
pothetical pairing was discussed[l[. After that, Vojta 
and Dagotto pointed the possible realization of odd- 
frequency spin-triplet s-wave (OTE pairing state) on a 
triangular lattice in the Hubbard model Recent de- 
tailed calculation by Yada [1] has supported this re- 
sult. Balatsky and Abrahams proposed an odd- frequency 
spin-singlet p-wave pairing (OSO pairing state) [4, l5|. 
Odd-frequency pairiiig has been studied on the Kondo 
lattice model [a, 0, Si- There are some experimental 
reports d which are consistent with the realization 
of the odd-frequency spin-singlet p-wave superconducting 
state (OSO pairing state) in Ce compounds [ol. [lOl [Tl| . 

It has been clarified recently that odd-frequency pair- 
ing correlation, i.e., pair amplitude is generated in inho- 
mogeneous system like superconducting junctions [T^fisl. 
H, IHj [H) 113 O'" vortex core [H, [ll| . It has been pointed 
out that OTE pair amplitude is induced in ferromagnet 
/ superconductor junctionsjlO, HH, [H, HI [H, HE [Hi and 
diffusive normal metal / spin-triplet odd-parity super- 



conductor iunctionsflSl. [la|. 

Stimulated by these preexisting works, we explore the 
situation, where odd-frequency pairing is stabilized. In 
this paper, we focus on the superconductivity on a quasi- 
one-dimensional triangular lattice. In order to clarify 
the dominant pairing state there, we solve the linearized 
Eliashberg's equation in the Hubbard model using the 
random phase approximation (R PA) and the fluctuation 
exchange (FLEX) approximation [11 HI HI]. 

We have found that odd-frequency spin-singlet p-wave 
becomes dominant pairing on a quasi-one-dimensional 
triangular lattice. In particular, this pairing becomes 
prominent on an isosceles triangular lattice due to geo- 
metrical frustration. 



II. MODEL AND FORMULATION 

We start with a single-band Hubbard model on an 
anisotropic triangular lattice as shown in Fig. [1] where 
tx, ty and ^2 are transfer integrals along x, y and diagonal 
directions, respectively. 




FIG. 1: A two-dimensional triangular lattice with transfer 
integrals t^, ty, and t2. 



The Hamiltonian is given by 
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where Uj denotes the transfer integral between two sites 
i and j. is the combination of nearest and second- 

nearest neighbors. (cia) and = c^q are creation 
(annihilation) and number operators, respectively. U is 
the on-site Coulomb repulsion. The band dispersion is 
given by 

Sk = —2tx cos kx — 2ty cos ky — 2t2 cos {kx + ky) . (2) 

In this study, the number of electrons per site is fixed to 
unity (half- filling) . 

In this paper, we calculate Green's function GiiSmk) 
in two different ways; using (i) the RPA and (ii) the 
FLEX approximation. 

(i) In the RPA, the Green's function is given by 
G(ie„, k) = {lEn — Efc + where £„ is the Matsubara 

frequency of fermion given by e„ — {2n — l)7rT with an 
integer n and fi is the chemical potential. Using Green's 
function, the irreducible susceptibility is obtained as 



T 

'n 



G(i£„, k)G{i£n + iuj,n, k + q), (3) 



where Um is the Matsubara frequency of boson given by 
ujjn = 2m'!TT with an integer m and N is the number of 
sites. The spin susceptibility is given by 



Xs(jt^m,g) = 



1 -UxoiiujTn,q)' 

and the charge susceptibility is given by 

Xc(«t^m,q) = -t—Tt — r \- 



(4) 



(5) 



(ii) In the FLEX approximation, first the bare Green's 
function Go{ien, k) = (je„ — efc -I- m)"^ is calculated. By 
substituting it into Green's function in Eq. ([3]), we ob- 
tain the irreducible susceptibility. The spin and charge 
susceptibilities are given by Eqs. ^ and (O, respec- 
tively. Using the susceptibilities, the effective interaction 
is given by 



3 1 

Vn{iuj,n, q) = -^U'^Xs{i(^r,i, ?) + -U'^Xc{i(^m, ?) " U'^Xo{it^m, ?)■ 



(6) 



r 



Then, we calculate the self-energy 



J:{ien,k) = —^Vn{iuJni,q)G{ien-iuj-m,k-q){7) 



N 



approximation, we solve the linearized Eliashberg's equa- 
tion. The effective pairing interactions for spin-singlet 
and spin-triplet channel are given by 



After using the Dyson equation 

G-l(^£„,fc) = Ga\ien,k)-^{isn,k) (8) 

we obtain new Green's function. The self-consistent it- 
erations are repeated until the sufficient convergence is 
attained. Thus, we obtain Green's function in which the 
self-energy is taken into account. 

Using Green's function, the spin susceptibility, and the 
charge susceptibility obtained in the RPA or the FLEX 
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V^{iu}^, q) ^ U + -[/^Xs(iw„, q) - -U'^Xc{iuJ,n, ql9) 

-^U^Xs{iuJ,n,q) - ic/^Xc(jWm,q),(10) 



V^{iujm,q) 



respectively. By substituting them into the linearized 
Eliashberg's equation for spin-singlet (spin-triplet) chan- 
nel 



AA(ze„,fc) = -jjYl Va''*Hi£n - iem,k - k')Giiem,k')G{-ie,n,-k')Aiie^,k'), 



(11) 



r 



the gap function A(ie„, k) and eigenvalue A are obtained. 
Superconducting transition temperature Tq corresponds 
to the temperature where A reaches unity. Thus, we con- 
sider that the larger the value of A becomes, the more 



stable the superconductivity becomes. In this paper, 
we take N = x Ny = 128 x 64 fc-point meshes. 
The Matsubara frequencies e„ and tUm have values from 
-{2Na - l)nT to {2Na - l)nT and from -2N^t:T to 
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2NcTtT, respectively, with Nc = 2048. 



magnetic spin fiuctuation[29'| since the superconductivity 
appears near the SDW phase. 



III. RESULT 

First, we study superconducting state on a quasi-one- 
dimensional triangular lattice based on the RPA. In our 
model, a quasi-one-dimensional triangular lattice is rep- 
resented by choosing the values of ty and t2 much smaller 
than that of t^- The quasi-one-dimensional direction is 
parallel to x-axis. 

In order to clarify how symmetry of gap function de- 
pends on the dimension of the lattice, we gradually 
change the lattice structure varying ty and t2 from a 
two-dimensional regular triangular one {t^ = ty) to a 
quasi-one-dimensional triangular one (t^ > ty). Here, 
we choose ty = t2 for simplicity. In this case, the lat- 
tice structure is an isosceles triangular one, in which spin 
frustration in the directions of ty and t2 exists. A of 
each pairing state changes as shown in Fig. [21 where 
Stoner factor {St = Ux^"'') is fixed to 0.97 by tuning 
the value of U. Here, Xo^"^ denotes the maximum value 
of xo- We see that OSO pairing state is dominant on 
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FIG. 2; (Color online) ty dependence of A for each pairing 
state, fixing ty = t2, obtained by the RPA with St = 0.97 and 
T/t^ = 0.04. 



a quasi-one-dimensional triangular lattice (ty/tx ~ 0.1), 
while ESE pairing state is dominant in almost all of the 
region. Hereafter, we focus on OSO pairing state. 

In order to study details of the OSO pairing state on 
a quasi-one-dimensional triangular lattice, we calculate 
the temperature dependence of A for ty/t^ — t2/tx = 0.1, 
and U/tx = 1.6. As shown in Fig. [3l A of OSO pairing 
state increases much above unity at low temperatures for 
U/tx ~ 1.6, where the Stoner factor St becomes almost 
unity. The value of St reaches unity at T/t^ ^ 0.05. This 
means that the OSO pairing state can be realized near 
the spin-density wave (SDW) phase. As in the case of 
ESE pairing, OSO pairing are are mediated by antiferro- 
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FIG. 3: (Color online) Temperature dependence of A for each 
pairing state by the RPA with ty/tx = t2/tx =0.1 and U/tx = 
1.6. 



The momentum dependences of the gap functions for 
the ESE and the OSO pairing states on a quasi-one- 
dimensional triangular lattice are shown in Fig. HJa) and 
(b), respectively, where the dashed and solid lines and 
the arrows represent the nodes of the gaps, the Fermi 
surfaces, and nesting vectors Q = (tt, 7r/2), respectively. 
We find that the momentum dependences of the ESE and 
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FIG. 4: (Color online) Momentum dependences of the gap 
functions for (a) ESE and (b) OSO pairing state by the RPA 
with ty/tx = t2/tx = 0.1, U/tx = 1.6 and T/tx = 0.06, 
where the dashed and solid lines and the arrows represent 
the nodes of the gaps, the Fermi surfaces, and nesting vectors 
Q = {n,n/2), respectively. 



the OSO gap functions can be approximated by cos 
(d-wave) and sinfe^; (p-wave), respectively. Here, we de- 
note d(p)-wave since this gap function changes sign four 
(two) times on the Fermi surface. It is noted that this 
p-wave has no nodes on the Fermi surface in the case of 
quasi-one-dimensional lattice {ty < 0.6). The shape of 
the Fermi surface becomes two lines with kx = ±7r/2 in 
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the absence of ty and t2- These two hnes are bent into 
an S-shape by introducing ty and ^2- It is noteworthy 
that these two "Fermi hnes" are perfectly nested with a 
vector Q — (7r,7r/2) at half-filhng with ty = t2- Because 
of this, the spin susceptibihty Xs{i^rm q) at q ~ (7r,7r/2) 
becomes strong as shown in Fig. El 




In the following, we study OSO pairing state focusing 
on frequency dependence. Matsubara frequency depen- 
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FIG. 6: (Color online) t2 dependence of A for each pairing 
state by the RPA with ty/t^ = 0.3, St = 0.995 and T/t^ = 
0.05. 
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FIG. 5: Momentum dependence of the spin susceptibility by 
the RPA with ty/t^ = t2/t^ = 0.1, U/t^ = 1.6, and T /t^ = 
0.06. 



It is known that on a two-dimensional triangular lat- 
tice, the antiferromagnetic fluctuation works equally in 
the directions of x- and y-axes. Therefore, in real space, 
neighboring two electrons in these directions with an- 
tiparallel spins make a Cooper pair. As a result, spin- 
singlet dx2-y2-'wave (ESE) becomes dominant. On the 
other hand, the antiferromagnetic fluctuation along x- 
axis becomes dominant on a quasi-one-dimensional tri- 
angular lattice. Then neighboring two electrons along 
cc-axis with antiparallel spins make a Cooper pair. As 
a result, spin-singlet d-wave (ESE) and p^-'wave (OSO) 
make pairs in x-direction. 

The gap function for ESE pairing state has nodes on 
the Fermi surfaces, while that for OSO pairing state is 
full-gap on the Fermi surfaces. This point is relevant 
to following fact that OSO pairing state dominates over 
ESE pairing state for sufficiently small magnitude of ty 
and ^2- 

We also explore the lattice with ty ^ t2 to clarify that 
ty = t2 is an indispensable condition for the realization 
of OSO pairing state. We gradually change the value of 
i2, fixing ty/tx = 0.3. The resulting A of each pairing 
state changes as shown in Fig. [5] We see that OSO pair- 
ing state is relatively-enhanced compared to ESE pairing 
state especially in the case of ty = ^2- This result is ro- 
bust against changing value of ty. The condition with 
ty = t2 corresponds to geometrical frustration, which 
suppresses the antiferromagnetic fluctuation in these di- 
rections. 



deuces of the gap function Aoso(*£n7 k) for OSO pairing 
state at fc = (''■/2, 0) and the effective pairing interaction 
V^{iWm,Q) with Q = (tt, 7r/2) for spin-singlet channel 
on a quasi-one-dimensional triangular lattice are shown 
in Fig. [7l Near the SDW phase, V^{iijJrm Q) has a sharp 




FIG. 7: (Color online) Matsubara frequency dependences of 
the gap function Aoso(ien,fe) at fe = (7r/2,0) and the ef- 
fective pairing interaction Va(iuLim,Q) at Q = (7r,7r/2) for 
spin-singlet channel by the RPA with ty/tx = t2/tx ~ 0.1, 
St = 0.95, and T/t^ = 0.05. 



peak at ujrn = in frequency space. After a simple trans- 
formation of the linearized Eliashberg's equation (fTT|) , we 
obtain following relation. 
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T T,n,m,k,k' Va^*\'i£n " iSm, k - k')G(iem,k')G{-iem, -k')A{iern, k')A{ien,k) 



InEq. lED, E„,fe|A(«en.fc)l'' G(^e™,fe')G(-^£™,-fc'), 
and (iSn — iSm, k — k') are always positive. Then nega- 
tive (positive) sign of Aoso(*^m, 's')Aoso(*£n: fe) niakes 
positive (negative) contribution to A. Due to the sharp 
peak of V^{iuj^a^Q) lj^^ = 0, pair scattering from e,„ 
to e„ with Em = £n makes main contribution to A. The 
gap function Ao^Q(iemk) changes sign in the process 
of scattering from k' to k through the nesting vector 
Q = (7r,7r/2), which makes the main contribution to A, 
in the momentum summation of the numerator. How- 
ever, scattering process for 7^ suppresses the value 
of A since gap functions with positive and negative sign 
of En have opposite signs each other in OSO pairing. 

Next we calculate the value of A with the FLEX ap- 
proximation in order to reveal how the above results is 
changed by the self-energy. As in the case of the RPA, 
to clarify how the superconducting state depends on the 
dimensionality of the lattice, we gradually change the lat- 
tice structure from a two-dimensional regular triangular 
one {tx — ty) into a quasi-one-dimensional triangular one 
{tx > ij/), fixing ty — t2. A of each pairing state changes 
as shown in Fig. [H Compared to the result of the RPA 
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FIG. 8: (Color online) ty dependence of A for each pairing 
state fixing ty = t2 by the FLEX approximation with 5** = 
0.97 and T /t^ = 0.02. 

in Fig. [21 the values of A are reduced. However, similar 
to the case of the RPA for ty/tx ^ 0.1, the OSO pairing 
state still remains dominant. 

We calculate the temperature dependence of A for 
ty/tx = t^jtx — 0.01. As shown in Fig. [9l the value of 
A for OSO pairing reaches up to 0.8. The present value 
of A is considerably high as compared to the correspond- 
ing values studied by the FLEX approximation in other 



I 

strongly correlated svstemsfsol [Sll [s^ [H, [11]. Up to 
now, only the values of A obtained for high-Tc cuprates 
and exotic systems with disconnected Fermi surface ex- 
ceed over the present value [H, [H, 113, HI, [H, SOj SH- 
Momentum and frequency dependences of gap function 
are qualitatively similar to the one obtained by the RPA. 
As regard the frequency dependence, there is only quan- 
titative difference. 
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FIG. 9: (Color online) Temperature dependence of A for 
each pairing state by the FLEX approximation with tyjtx, — 
t2/t^ = 0.01 and U/t^ = 2.5. 



Matsubara frequency dependences of value of 
G{ien,k)G{—ien, —k) at k = {tt/2,0) obtained by the 
RPA and the FLEX approximation, which directly affect 
the value of A as noted in Eq. (fT2|) . arc shown in Fig. 
[TOl G{ien, k)G{—ien, —k) for the FLEX approximation 
is smaller than that for the RPA since imaginary part 
of self-energy corresponds to damping of quasi-particles. 
Thus, the presence of self-energy decreases the value of 
A in the FLEX approximation. 

Matsubara frequency dependences of the normalized 
effective pairing interactions V^{iuJm,Q)/V^{iO,Q) at 
Q = (tt, 7r/2) for spin-singlet channel in the RPA and 
the FLEX approximation are shown in Fig. [TlJ A 
peak width of V^{iujm,Q) for the FLEX approximation 
is broader than that for the RPA in Matsubara frequency 
space due to the presence of the self-energy. This broad- 
ness of the peak width of the effective pairing interaction 
V^{iujm, Q) relatively enhances the value of V^{iuJrmQ) 
with LOm 0. Subsequently, the scattering processes 
from positive e„ to negative e„ relatively increase the 
summation of Matsubara frequency in the numerator of 
Eq. (112). Thus, the OSO pairing state is suppressed. On 
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FIG. 10: (Color online) Matsubara frequency dependences of 
G{i£„,k)G{-ie„,-k) at fe = (vr/2,0) by the RPA and the 
FLEX approximation with ty/t^ — t2/tx = 0.1, St ~ 0.95, 

and T/t, 
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FIG. 11: (Color online) Matsubara frequency depen- 
dences of the normalized effective pairing interactions 
Va{iujm,Q)/Va{iO,Q) at Q = (7r,7r/2) for spin-singlet chan- 
nel by the RPA and the FLEX approximation with ty/tx = 
t2/U = 0.1, St = 0.95, and T/U = 0.05. 



by solving the linearized Eliashberg's equation based on 
the RPA and the FLEX approximation. Surprisingly, 
odd-frequency spin-singlet Pa:-wave (OSO pairing state), 
which is not familiar, is enhanced near the SDW phase. 
The OSO pairing state becomes prominent on an isosceles 
triangular lattice {ty = t2), where the geometrical frus- 
tration is the most significant. Even if the self-energy is 
introduced by the FLEX approximation, above conclu- 
sions are not changed. 



The OSO pairing state is induced in the following way. 
In real space, neighboring two electrons with antiparallel 
spins make a Cooper pair mediated by the antiferromag- 
netic spin fluctuation near the SDW phase. The value 
of the effective pairing interaction Vi^{iiUm,Q) for spin- 
singlet channel at the nesting vector Q = (tt, 7r/2) has a 
sharp positive peak at uim = in the frequency space. 
In this case, pair scattering with conserving frequency 
makes major contribution to the value of A. Therefore, 
sign inversion of gap function through the nesting vector 
Q — (tt, tt/2) in the momentum space enhances the value 
of A in Eq. (fT2|) because the effective pairing interaction 
for spin-singlet channel has a positive value. This fa- 
vors Pa;-wave, which is full-gap on the Fermi surface. In 
accordance with Fermi-Dirac statistics, spin-singlet p^- 
wave can be interpreted as odd-frequency pairing. These 
results in this paper suggest possibility of odd-frequency 
superconductivity realizing in bulk system. 



In the present paper, only on-site Coulomb interaction 
is considered. In the presence of off-site Coulomb interac- 
tion, it is known that SDW and CDW can compete each 
other in quasi-one-dimensional superconductor. In that 
case, competition between even-frequency spin-triplet /- 
wave pair and even-frequency spin-singlet rf-wave pair has 
been pointed by several theories [H, |43, 0, |40| . It is 
a challenging issue to consider the possible realization of 
odd-frequency pairing in the presence of off-site Coulomb 
interaction. 



the other hand, for the ESE pairing state, the above scat- 
tering processes do not suppress the value of A since the 
gap function has always same sign in frequency space. As 
a result, the OSO pairing is suppressed more significantly 
by the self-energy as compared to the ESE pairing. 

IV. CONCLUSION 

We have studied symmetry of gap functions on a quasi- 
one-dimensional triangular lattice in the Hubbard model 



Beside this problem, to clarify the superconducting 
properties of odd-frequency superconductor is an im- 
portant problem. Since phase sensitive probes, e.g., 
tunneling and Josephson effects, are crucial to identify 
the pair in g symmetry in unconventional superconductors 
[13,111, EmI' similar studies on odd- frequency supercon- 
ductors become important [H, \E^, HH, . It is necessary 
to calculate temperature dependence of energy gap func- 
tion to reveal the superconducting properties. 
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